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ABSTRACT 


Characteristic matrix method has been used in this work to determine the mode field 
profile and its evolution in multilayer planar waveguides. The fields in any one layer is 
expressed as multiplication of the characteristic matrix to the field at its neighboring 
layer, ultimately relating it to the fields in the substrate layer. Two symmetrical parallel 
optical waveguide couplers considered that are separated by a center-to-center distance 
s. The power transfer distance under which the power is completely transferred from 
wave-guide 1 to wave-guide 2 called coupling length ‘Lc’ has been calculated using the 
normal modes of this multilayer structure. Dependence of ‘Lc’ with the thickness of the 
film (0 rib hight (h) rib width (w), waveguide spacing ( 5 ) and refractive indices if 
substrate (n^) and the film (r/) have been studied by this technique. Tapered optical 
waveguide couples have also been studied by this technique. In the first case the input 
rib waveguide was tapered (0) in shape so that width of this waveguide monotonically 
reduces with distance, whereas the coupled waveguide was kept normal. The separation 
between the waveguides were kept constant along the length of the coupler. In the 
second case the input tapered waveguide reduces in width along the length of the 
coupler, whereas the coupled waveguide increases in width along the length. Single 
mode planar waveguide structure were thus design for a rib waveguide use is made of 
the effective index method to formulate a multilayer planer waveguide and the 
characteristic matrix method is repeated to obtain the characteristics of the lowest order 
mode of the rib waveguide. The three cases of coupler i.e. (i) parallel coupler (ii) One 
waveguide tapered coupler, and (iii) Both waveguides tapered coupler were coupled 
with respect to its coupling lengths and insertion losses. The coupling lengths were 
lower for the single tapered case with respect to that of the normal coupler, and for the 
both waveguide tapered was lower compared to the single waveguide taper. For the both 
waveguide taper coupler the coupling length decreased with the taper angle (0) and also 
with the separation (.y). The separation distance was taken to be 5 < w . However, 
separation had little effect on the coupling length where 0 was beyond O.SMnsertion 
losses in the three cases are 0.1873 dB, 0.9404 dB and 6.4588 dB respectively. 
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Chapter 1 


INTRODUCTION 


Film waveguides are the basic structure for both the passive and the active devices in 
integrated optical circuits. The physical principles involved in these waveguides and 
other related thin-film structures, coupled with those of modem laser electronics, form 
the foundation of integrated optics. An optical film waveguide is simply a sheet of thin 
film deposited on a substrate, which has a refractive index, larger then substrate 


refractive index _ the air-space or cladding refractive index is Such a three- 
layered arrangement, with the highest refractive index in the middle layer, forms a lens 
like medium, which has a natural tendency of focusing a light beam into the film. 
Because of this focusing action, it is not surprising that a thin film can be used as a 
waveguide. Consider fig. 1.1 where a ray of light starting in the film and propagating 
toward the upper surface of the film is total internally reflected at the interface of the 
film and air and the film and substrate. This structure, known as slab waveguide, guides 
the light beam through the film. 
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Figure 1.1 Wave propagation in slab waveguide. 


Analysis of the field profile of a guided wave is important to understand the coupling of 
a lightwave from its source to a waveguide and from one waveguide to another 
waveguide. The mode pattern is the standing wave pattern produced by the interferences 
in the transverse direction of the optical guide. Some waveguides are designed to have 
certain mode profile to provide particular spot shape and better coupling efficiency 



between different types of waveguides. Moreover, measuring the mode profile is one of 
the popular methods for characterizing waveguides in practice. Analysis of the mode 
profiles of a 3-layer slab waveguide, which provides, such guide characteristics as the 
effective guide index and effective guide width, is dependent on waveguide parameters. 
These are the refractive indices and ric of the substrate, the film and the cover 
respectively. The film thickness w and the free space wavelength X or propagation 
constant, ko=2Tz/X =co/c of the light. The number of independent parameters can be 
combined by the introduction of appropriate normalized parameters to provide universal 
plots from which the effective index and effective guided thickness can be determined 
for any slab waveguide configuration, this will be discussed in chapter 2. For basic 
guide parameter the normalized frequency or normalized film thickness V, defined as: 

V = k^w^n] -n] 

The mode profile of a 3 -layer slab waveguide can be carried out in a straightforward 
maimer [1], special methods have been developed for more complex structures. Well 
known among these are the characteristic matrix method [2], finite element method, 
Fourier expansion method, and discrete sine method. 

X 



Figure 1 .2 Sketch of a multilayer slab waveguide with substrate index and cover 
index ric. They-axis indicates the direction of mode propagation. 



The characteristic matrix method has been used to analyze the coupling in thin film as 
well as the analysis of multilayer planar waveguides. Here the characteristic matrix 
method is used for its application to study the mode profile as well as its evolution in>'- 
invariant planar waveguides with arbitrary refractive index profile. This method will be 
shown to be relatively simple and requires only modest computational power. Guided 
mode fields extend outside the guiding region into the surrounding cladding medium, 
called evanescent field, decaying exponentially away from the waveguide cladding 
boimdary. For two waveguides close together such that the evanescent fields of the 
propagating modes in two waveguides overlap, there exists coupling of optical power 
between the two waveguides. Such a device is called a directional coupler. The fraction 
of power coupled per unit length from one waveguide to the other is determined by the 
overlap of the modes between the two waveguides as shown in fig. 1.3. 


X 



Figure 1 .3 Electric field profile in symmetrical coupler. 


Thus it depends on the separation of the guides, the mode penetration into the substrate, 
and the interaction length. This energy exchange can be used in building an optical 
modulator or optical switch. Fig. 1.4 shows a directional coupler consisting of two single 
mode waveguides, which are in close proximity over a length Lc. When the incident 
light is focused into a single channel at y=0, then there is a periodic exchange of energy 
between the two waveguides. If the propagation constants of the modes in the 
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waveguides are identical then there is a complete transfer of energy from one 
waveguide to the other over a length “Lc” and if the propagation constant is not identical 
then there is only an incomplete transfer of energy for the same length i.e. the coupling 
length Lc become different. 

y=0 



Figure 1 .4 Optical directional coupler in which two waveguides are at close proximity 

over a length Lc- 


A theoretical analysis of the coupling between two dielectric waveguide usually 
involves a formal solution of Maxwell’s equation using matrix method, when the 
waveguide modes have significant overlap. The waveguides can have any arbitrary 
cross-section however for the sake of simplicity only guides with rectangular cross- 
sections have been considered here. In the following chapters, chapter 2 discuss the 
matrix method used to establish the mode field profiles and propagation constants of 
multilayered slab waveguides. Chapter 3 discusses the analysis of multilayer rib 
waveguides using the effective index method. Hence the quasistatic modal field and 
propagation constant for a tapered waveguide along its length have also been discussed. 
Coupled mode analysis have been carried out in chapter 4 for parallel waveguide 
couplers, however a matrix method is developed for asymmetric as well as tapered 
waveguide coupler using the lowest order mode in the effective index method. The 
result of all the analysis is presented in chapter 5. Hence the field profiles as well as 
power profiles of tapered waveguide couplers are presented. From this the coupling 
lengths of tapered couplers are obtained and compared with parallel waveguide 
couplers. Chapter 6 concludes this work. 
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Chapter 2 

Multilayer Slab Waveguide 


2. 1 Characteristic Matrix Method 


A z-invariant multilayer planar waveguide with p+1 layers is shown in fig. 2.1. This 
structure can represent either step index waveguides or a stepwise approximation to 
graded index waveguides. The waveguide medium is assumed to be linear, lossless, 
source free and non-magnetic. Light is propagated along positive z direction. In order to 
simplify the analysis, TE and TM polarization field will be considered separately, in 
rectangular coordinates system field vector F(E or H) expressed as 

F = Fix,y)e-^^ (2.1) 


where P is called propagation constant of guide 
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Figure 2.1 Multilayer planer waveguide structure 


2.1.1 Transverse and Axial fields 


The first two Maxwell’s equations are represented as 


V X if = (cr + jcos)E 

(2.2) 

V xE = -jco/uH 

(2.3) 


where o and e represents the conductivity and permittivity of the dielectric respectively. 



0 


For a rectangular coordinate system with o =0, Maxwell’s Eq. (2.2) and Eq. (2.3) yield 
the scalar equations 


+ 

II 

H 

1 

1 

(2.4) 

ax 

(2.5) 

dx dy 

(2.6) 

+ 

II 

(2.7) 


(2.8) 

dx dy 

(2.9) 

2 . 1 .2 TE Polarization Mode 


For TE polarization, Ex=0, Ez=0 and Hy=0, 

the Maxwell’s equations lead 

relations 


coil 

(2.10) 

II 

1 

1 

N 

(2.11) 

- + j/3H^ = -js^aE 

OX 

(2.12) 


By defining new variables U = mAY = co/uH^ , the equations are converted to 


(2.13) 

(2.14) 
itg = 2n:/Ji 


U' = -jY 
U" + k^U = 0 

where U' = 


; 



and n is the refractive index of the layer under consideration 
Solutions of Eq. (2.13) and Eq. (2.14) in the layer have the general forms 


/ 


U{x) = A tx'p{-jkx) + B exp(yfcc) (2.15) 

Y {x) = k{A exp(-jrfoc) - B exp(yfoc)] (2.16) 

By applying boundary conditions, U and Y at any interface plane can be related with its 
adjacent planes by means of the characteristic matrix associated with the layer in 
between expressed by 


U. 

Y 


= Mi 






(2.17) 


where the characteristic matrix M, is given as 
M. 


cos{k^di ) - j sin(A:,d, ) / A, 

-jK^sinik^,) cos{k,d,) 


(2.18) 


Eq. (2.17) is generalized to give the relation 


r^,i 




= M, 



7o. 


where Mr - MiMi.i M 2 M 1 = 


m, 


rri: 


m, 


(2.19) 


( 2 . 20 ) 


If decay constant in substrate is ys then, Eq. (2. 1 5) and Eq. (2. 1 6) can be written as 
t/,(x) = C/oexp(x,x) (2.21) 

Y.Xx) = jy.U. ( 2 - 22 ) 

where - n]kl at interface 0, we obtain 

C/o=f/.v(0) (2.23) 

Yo=VM = jrsUo (2-24) 


In terms of Uo and Yo, the values of 17 and F at any point x, between x, and xt+i can be 
obtained by envisioning a virtual interface at x, according to 



'Uix) 


"cos(^,,, (x - X, )) - ,/sin(A,^,(x - 

/(x)_ 


- sin(x-,„ (x - x, )) cos(A:;^, (x - x,. )) 


X 


^0 


which result in U (x)^ 

where 


(2.25) 

(2.26) 


= ("J,„ + )cos(A:,,,(x - jc,)) - 


y(^,„ +J’rs^i ,, ) sin(^,^, (;c - x, )) 


\+i 


(2.27) 


If decay constant in cover or cladding is yc then the value of U and Y in the cover can 
similarly obtained by writing 

t/,(x) = 5,exp(-y,x) (2.28) 

lXx) = -jnU, (2.29) 

where y] = - nlk] at the interface p, we have 

rr/.i 

(2.30) 


^p' 


'5,exp(-7^,x^) 


'u; 

Yn. 


-jyA^^v{-rc^p)_ 

p 

/o. 


which result in 


U,{x) = gexp(-y,,x)f7o 


(2.31) 


where 


g = 




Vrc^wi-TcXp) 

From Eq. (2.30), we obtain the dispersion relation for TE modes 
jX/s^p,, + r.^p „ ) = 


(2.32) 


(2.33) 


For TE modes parameter “b ” is called the normalized guide index [1] and which is 
defined by 



{n]-n]) 


( 2 . 34 ) 


where n^jf is the effective index of waveguide, defined as jS/ko. 

2. 1 .3 TM Polarization Mode 

For TM-polarized mode, Hx-Q, Hz=Q and £y=0, the following notation is introduced 


U = H.. 


( 2 . 35 ) 


Y = coSqE, 


( 2 . 36 ) 


This leads to a solution of the wave equation as 


U{x) = Aexp(- jkx) + Be.xp{jkx) 


( 2 . 37 ) 


T(x) = — Y{Atxp{-jkx) - 5exp(7foc)] 
n 


( 2 . 38 ) 


The characteristic matrix is then given by 


cosijk^d^ 

- jk^sinikid^) I nf 


-jnfsm(k,d,)lk, 

cos(k,di) 


( 2 . 39 ) 


For the mode fields in substrate again 


^.v(^) = ^0 exp(r,x) 


( 2 . 40 ) 


=~-TTvt^oexp(r.vX) 


( 2 . 41 ) 


Following the same formulation described for the TE case, one may write for TM 




= f(x)x,^xsx,^,^0 


( 2 . 42 ) 






( 2 . 43 ) 
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Similarly for the field in cover layer, we have 

t/,(x) = gexp(-y^x)f/o (2.44) 

with 



For a waveguide with known geometric structures, excited at a known wavelength, the 
propagation constant P and the number of guided modes for TE and TM polarization 
may be obtained by solving Eq. (2.33) and Eq. (2.46) respectively. Hence, the electric 
field profile Ey{x) in substrate, cladding and central layer for TE polarization mode can 
be obtained from Eq. (2.21), Eq. (2.31) and Eq. (2.26) respectively while the magnetic 
field profile F/^(x) in substrate, cladding and central layer for the TM-polarization mode 
can be obtained from Eq. (2.40), Eq. (2.44), and Eq. (2.42) respectively. 



Chapter 3 

Waveguide Analysis 


3.1 Rib W aveguide 

For a N-layer rib waveguide, the refractive index values (substrate), m, ... on 
(L ayers 1 to N), (cover), the thicknesses ti,..., of the inner layers, the etching 
depth h, and the width w of the rib. All dimensions are meant in micrometers. The 
figure illustrates the relevant geometry 
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Figure 3.1 Multilayer rib waveguide structure 


Note the orientation of the transverse coordinate axes with the x-direction perpendicular 
to the film plane. Light propagates along the z-direction. The entire structure is meant to 
be symmetric with respect to the x-z plane. The structure informs about propagation 
constants ^ (pm’') and effective mode indices rieff. The effective index simulation starts 
with computing the one-dimensional modes of the central slice -w/2 < z < w/2 and the 
two identical lateral slices z < -■w/2 and w/2 < z. 


3.1.1 Effective Index Method 

In this method, the effective index (El) of the structure is obtained by successively solving 
two transcendental slab equations. If the example of a rib waveguide is considered, fig.3.2 
(a), then the method in the first step solves the transcendental equations for three vertical 


slabs as shown in fig.3.2 (b). The effective indices Ueffi and so obtained are used as 
the refractive indices for a 3-layer horizontal slab waveguide, as shown in fig.3.2(c). 
Solving the transcendental equation for the horizontal slab waveguide gives a good 
approximation to the effective index of the original rib waveguide structure, for the lowest 
order mode. The advantage of the El method is that it can be applied to a wide variety of 
structures. 



(c) 

Figure 3.2 The Effective Index method for the rib waveguide (a) the original rib wave 
guide, (b) solving the vertical slab problem to define We#; and (step 1), (c) 
equivalent horizontal slab structure (step2) 

For TE modes the normalized guided indices b is defined by 




(3.1) 


where n^ffi . neff 2 and rieffSXQ the effective index of inner, outer regions and overall structure 
respectively. 

The normalized frequency or film thickness V for final structure will be 

V=ko w ( 3 . 2 ) 


where w is the width of central region and ko=2n/X. 


The mode profile plots show the field of the basic electric component Ey for TE-like 
modes and of the basic magnetic component Hy for TM-like modes. Note that the 
effective index approximation assumes factorizing mode fields one due to lateral and 
another are due to central slice. Hence the horizontal profile yields the correct y- 
dependence for all positions x. According to the initial decomposition in vertical slices, 
there are two different vertical profiles, one corresponding to the central region and one 
for the lateral slices. These curves are the constituting 1-D mode profiles. 

3.2 Quantum Well Waveguide 


When dimension of the semiconductor becomes comparable to the mean free path (mfp) 
of carriers, quantum size effects become important and they dominate the electronic 
properties of the material. Of particular interest are quantum well stmctures [3], show in 



Figure 3.3 Schematic representation of a quantum well in three dimensions formed by 
two semiconductors with bands Egl (InGaAs) and Eg 2 (GaAs). d] and d 2 are well and 

barrier widths respectively. 


In such structures two semiconductors of different band gap energies alternate to form 
a synthetically modulated structure. If the narrow gap semiconductor layer thickness 
is comparable to the m^ of electrons in the material, the motion of carriers in the 
direction perpendicular to the interfaces is quantized. If the thicknesses of the wide- 
bandgap layers are thin enough, such that carriers may tunnel through, the discrete 



energy levels formed in the quantum wells due to size quantization splits into 
minibands. Such a composite material is called a superlattice, which has electronic 
and optical properties that are unique to the whole structure. For example, carriers are 
transported through such material by miniband conduction. Also, the effective 
bandgap of entire composite structure is approximately equal to the energy separation 


A InxGaasi-x =25 A, GaAs=70 A; 

♦ InxGaasi.x=50 A, GaAs=100A; 

T InxGaas i .x = 1 00 A, GaAs=5 0 A; 

■ lUxGaasi-x =70 A, GaAs=25 A; 

# InxGaas i-x =30 A, GaAs=45 A. 

The insert show an enlarge view of the 
characteristics around x =0.2. 


Figure 3.4 Refractive index step between InxGai-xAs and GaAs for different values of 
X as indicated by the representative symbol [4] 

Usually in compound semiconductor materials with =\eV , such coupling between 

wells occurs for barrier thickness < 5 nm. For barrier thickness larger than 5nm, the 
individual well are largely uncoupled and the electronic properties in a single-quantum 
well (SQW) or multiquantum well (MQW) are quite similar. In general the refractive 
index of a strained InxGai-xAs is given as [4] 

ni=n2+An (3.4) 


between electron and hole miniband 



where An is shown in fig.3.4. 



So the refractive index a MQW of GaAs/InxGai.xAs is 

_ 


«/=■ 


+ d2 


(3.3) 


The value of n, for a specific structure as shown in fig.3.5 can be calculated using 
the dn [4]. 
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Figure 3.5 Schematic representations of n^GaAs quantum well structure. 


3.3 Tapered Waveguide 

In the given structure of a tapered waveguide fig.3 .6, if signal is incident at the end B 
(W 2 ) it will propagate through the structure and finally it reaches at the end A (w/) but 



Figure 3.6 Schematic representation of multilayer tapered waveguide. 


due to tapered structure of the rib its modal width changes and so does the 
propagation constant. The effective index method be used in the equivalent waveguide 



structure as shown in fig.3.7, where at different y the modal profile and the 
propagation constant will be different, under the quasistatic approximation. 



Figure 3.7 The equivalent waveguide representation of multilayer tapered waveguide, 
rieffs and rieffi are effective indices of lateral and central slice of the structure. 



Chapter 4 


Couplers 

4. 1 Normal Mode Analysis 

Instead ol coupled waveguide analysis, another approach of great practical interest is 
coupling of the normal in the steady state of the composite structure. For example, two 
parallel optical waveguides are coupled to each other via their evanescent fields. A 
wave set up initially in one guide is transferred to the other guide. Another coupling is 
affected between forward and backward waves by periodic perturbations on an optical 
waveguide. This process is also akin to the coupling of the normal mode analysis [5] 
presented here. 

z z z z 



Figure 4. 1 Slab optical coupler with symmetry in the refractive index distribution. 

Fig.4. 1 [6], correspond to the two lowest order modes of the pair of waveguides in a 
directional coupler considered as a single waveguide stmcture. Energy incident on one 
waveguide can be considered to be a linear combination of the two modes, which is 
shown in the above figure. As the modes propagate in the waveguides, they develop a 
phase difference due to difference in ^ values of the two modes. Thus at some point a 



phase difference of n is developed which corresponds to the total energy in the other 
waveguide. 

4.1.1 Basic formalism 

Consider two waves a/(y) and <t2(y), of modes 1 and 2, which in the absence of coupling 
have propagation constants /?/ and p 2 - They obey the equations 

-7'- = -'MW (4.1) 

dy 

da^ ^ . 

-r = -'M2(y) (4.2) 

dy 

Suppose next that the two waves are weakly coupled by some means, so that a\(y) is 
affected by a^iy) and a 2 (y) is affected by ai(y). Then, the equations become 


^ = -WMy) + Ky^a^iy) 
dy 

(4.3) 

dy 

(4.4) 


where the constants K 12 and K 21 represents the strength of interaction between the two 
modes and referred to as coupling constants, in the absence of any interaction between 
the two modes ki 2 = K2i=0. The above equations show that in the presence of interaction, 
the amplitude of the mode in a waveguide depends on the amplitude of the mode in the 
other waveguide. The coupling constants depend upon the waveguide parameters, the 
separation between the waveguides and the wavelength of operation. For identical 
waveguides one can show that 

^■21 = K]2 (4.5) 

(4.6) 


In order to solve this set of equations, the solution will be 
a,(j) = fl,(0)exp(-/^) 

«2(T) = «2(0)exp(-/^) 


(4.7) 



the existence of a wave in the system consisting of the two coupled waveguide, which 
propagates with a phase constant p and which is a superposition of the modes of 
waveguides 1 and 2 with amplitude a/(0) and a2(0). Substituting from Eq. (4.6) and 
(4.7) in Eq. (4.3) and (4.4) 

-//X3,(0) = -ip^a^{Q)-iK^^a^{Q) zXy=0 


= 0 (4.8) 

= 0 (4.9) 

solve the above two equations Eq. (4.8) and Eq. (4.9) 

- A A +/?,) + (^,/?, ) = 0 (4.10) 

where k = (x‘i 2 't' 2 i)''^^ ; iJ^us 

A..=^(A-A)±[.j(A-A)+'f’l''' («i) 


Thus, in the coupled waveguide, one has two independent sets of modes, one 
propagating with propagation constant Ps and the other with Pa. Therefore the general 
solution of Eq. (4.3) and Eq. (4.4) can be written as 

a, (y) = a, exp(-//?,y) + a„ exp(-//?„y) (4. 1 2) 


(y) = — a, exp(-Zy0,y) + — - - - a, exp(-z^„y ) (4. 1 3) 

K’i2 X',2 

It is noted that when the two waveguides are identical then Pi= P 2 =Po 


Ps=Po+K 

(4.14) 

Pa=P(rK 

(4.15) 

and we have for the mode with propagation constant P= Ps 


«7(0)=a2(0) 

(4.16) 


and for the mode with propagation constant P= Pa 


aiiQ))= -a2{Q)) 


(4.17) 



We now assume that at y=0, the mode in waveguide 1 is launched with unit power and 
that there is no power in waveguide 2.Then we have they=0 initial conditions. 


+a.=l 


A-/?,.. , A,-/?, 


which gives 


o, + 


A -A 
A - A 




' a =0 


and a„ 


(4.18) 

( 4 . 19 ) 


A - A- 
A - A 


The power in waveguide 1 and 2 is proportional to |a,(>')|^and |a 2 (>^)f respectively. 


«.(T)f =1- 


K(T)f = 


K 


{A/3l2y+K^ 


sin'[((A/?/2)^+/c')''^;;] (4.20) 


/c‘ 


(AjS/iy+K^ 


sin'[((AA/2)^+/rA''V] (4-21) 


where A/3 = A ~ A 2 


As can be seen from Eq. (4.20) and Eq. (4.21), there is a periodic exchange of energy 
between the two waveguides with a period 


/ = 


71 


{{Afiiiy+Kn 


2 \\n 


(4.22) 


Hence the length 


L = l/2 = 


It 

2{{AI3l2f+K^f^ 


(4.23) 


This is the minimum interaction length under which max power can transfer it is called 
coupling length of the directional coupler. From Eq. (4.23), it is clear that complete 
power transfer can take place when A y?=0, i.e. the propagation constants of the mode in 
each waveguide are same. Then coupling length will be 



K 


71 


(4.24) 
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Figure 4.2 Schematic representation of multilayer symmetric coupler. 



Figure 4.3 The equivalent waveguide representation of multilayer symmetrical 

coupler. 

4.2 Tapered Coupler 
4.2. 1 Single Tapered Coupler 


A typical corrugated single tapered waveguide coupler is shown in fig 4.4.This structure 
has two rib waveguides A and B separated by distance s, A is tapered rib waveguide of 
tapering angle 9 and waveguide B a parallel rib waveguide, the etch depth and refractive 
index of each of them are h and rip respectively, the initial (y=0) widths of waveguide A 
and B are wi and Wz Let refractive indices and film thicknesses of another p-l layers are 

njn 2 rip.] and tit 2 tp.\ respectively. Due to asymmetry in the guided structure A, 

its width and propagation constant changes with length y (y=0 to Lc), while waveguide 
B has constant value of width and propagation constant throughout the length. Effective 
index method can be applied to calculate the value of propagation constant and effective 


index ot guided structure at each point (y=0 to Lc), in the fig.4.5 the equivalent structure 
is shown. Assuming that a field is launched in waveguide A, which after the coupling 
length will be transferred to waveguide B. The normal mode fields are evaluated for the 
tapered waveguide structure in a similar manner to that of the previous section and 
coupling length are found. 
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Figure 4.4 Schematic representation of multilayer asymmetrical coupler in which one 
waveguide is parallel while another is tapered. 
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Figure 4.5 Equivalent waveguide representation of multilayer asymmetrical coupler 
4.2.2 Double Tapered Coupler 

A typical geometrical structure double tapered waveguide coupler is shown in fig 4.6. 
The structure is similar to single tapered coupler except both the rib waveguides A and 
B are tapered with tapering angle 0, the separation of waveguides is s, which will 
remain constant throughout the length. Again effective index (El) and characteristic 




matrix method, which is discussed earlier used to calculate coupling length of the 
waveguide, so the equivalent structure of the waveguide is converted into fig.4 


Interface p 


Interface 1 



Cover nr 


Substrate 


Figure 4.6 Schematic representation of multilayer asymmetrical coupler in which, both 

waveguides are tapered. 



Figure 4.7 Equivalent waveguide representation of multilayer double tapered coupler. 



Chapter 5 

RESULTS AND DISCUSSION 

5.1 Slab Waveguide Characteristics 

By the characteristic matrix method section.2.1, the complete dispersion relation in 
multilayer slab waveguides, the Eq. (2.33) and Eq. (2.46) represents for TE and TM 
guided modes respectively have been evaluated for a normalized propagation constant 
‘ P. for difterent guided structures the normalized wavelength dependent characteristics 
for TE and TM modes {h-V") curve are shown in fig.5.1, fig.5.2 and fig.5.3 for 4,5 and 6 
layers respectively. 



Figure 5.1 Plot of normalized propagation constant {b) against normalized frequency (V) 

for 4-layer slab waveguide. 
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2 Plot of normalized propagation constant (b) against normalized frequency (V) 
for 5-layer slab waveguide. 
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Figure 5.3 Plot of normalized propagation constant (b) against normalized frequency (F) 

for 6-layer slab waveguide. 





5.2 Electric field (£) vs Distance (x) 

The transverse electric (TE) and transverse magnetic (TM) field profile in the same 4,5 
and 6 layers slab waveguide structure are shown in fig.5.4, fig.5.5 and fig.5.6 
respectively. In all cases n=3.65 layer is assumed to be the guiding layer where the peak 
of the electric field lies. As expected, it is also observed that the evanescent tail is larger 
in the substrate than in the air. 
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Figure 5.4 Plot of electric field (£) against distance (x) fgr 4-layers slab waveguide. 
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5.5 Plot of electric field (E) against distance (x) for 5-layers Slab waveguide 
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Figure 5.6 Plot of electrical field {E) against distance (x) for 6-layers slab waveguide 






5.3 Single Mode Rib Waveguide Characteristic 

The effective index method, which is discussed in section 3. 1.1, we can calculate the 
effective index of central slice and lateral slice of rib structure now the whole structure 
is converted into three-layer slab waveguide, the effective indices of equivalent rib 
structure is shown in fig.3.2. 



Figure 5.7 Electric field (£) profile in rib waveguide. Solid line represents profile in 
central slice where dotted profile line represents profile in lateral slice. 



Electric field 'E' (norm unit) 


Figure 5.8 Electric field (E) distribution in the effective guided structure. 





The fig.5.7 and fig.5.8, which is given above represents the electric field profile of the 
real structure and equivalent structure of rib waveguide respectively. Once we get the 
effective guided structure we draw the normalized propagation constant (b) vs 
normalized frequency (V) characteristic, as shown fig. 5. 9, by this curve we can calculate 
the effective index of equivalent rib structure. 



Figure 5.9 Plot of normalized frequency (F) against normalized propagation 
constant (b) for effective guided structure. 

After calculating the effective indices of rib structure now we discuss the normalized 
propagation constant (6) vs width (w) of rib waveguide, fig 5.10 shows the profile for 
the given structure, by this figure it is clear that with increase of the value of width, 
normalized propagation constant increases. Another fig. 5.11 represents normalization 
propagation constant (b) vs wavelength (1) characteristic, it is clear from this 
characteristic that with increase of the value of wavelength normalization propagation 
constant (b) decreases. Fig.5.12 shows the normalized propagation constant vs etching 
depth (h) characteristic for the given structure, from this figure it is clear that with 
increase of the value of etching depth, normalized propagation constant (b) increases, if 




etching is keep on increasing a point come when lateral slice completely removed and 
characteristic become straight line. 



Figure 5.10 Plot of normalized propagation constant {b) against width (w) for single 

layered rib waveguide. 



Figure 5.11 Plot of normalized propagation constant (&) against wavelength (1) for 

effective guided structure 






Figure 5.12 Plot of normalized propagation constant (b) against etch height (h) for 

single layered rib waveguide 


The fig.5.13 shows the electric field profile in the rib structure, aty=0 signal is launched 
in the wave guide and it propagates throughout the length L, 



Figure 5.13 Electric field (£) distribution in single layered rib waveguide with respect to 

length (L) and width (w). 



it is noted that field profile does not changes throughout the length it is because of, the 
propagation constant ^ does not change with change in length. In fig.5.14 power 
distribution characteristic of rib waveguide is shown, it is clear from this figure that 
power distribution throughout the whole structure will remain constant, it will not 
varied with change in distance. 
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Figure 5.14 Distribution of power (P) in single layered rib waveguide with respect to 

length (L) and width (w). 

5 .4 Quantum W ell W aveguide 

The structure of quantum well is shown in section 3 .2, by the help of Eq. (3 .3) and Eq. 
(3.4) we can calculate the refractive index of the quantum well structure so the real 
structure converts into multilayer slab waveguide. The quantum well structure in the 
equivalent waveguide is work as a film and most of the signal propagates through this 
only. The figure, which is given below, shows the normalized propagation constant (b) 
vs normalized frequency {V) characteristic for equivalent quantum well waveguide, both 
he TE and TM propagating modes are represented by solid and dotted lines respectively. 
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Figure 5.15 Plot of normalized propagation constant (b) against normalized 
frequency (V) for quantum well waveguide. 
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Figure 5.16 Electrical field (£) distribution in quantum well waveguide at different 

distances (x). 






From the fig.5.15 it is clear that for TE and TM modes only two, zero order and first 
order mode will propagate in the given waveguide structure, the electric field in 
substrate, cladding and central layers can be calculated by Eq. (2.21), Eq. (2.31) and Eq. 
(2.26) respectively. The electric field profile for equivalent structure is shown in 
fig.5.16. 

5.5 Tapered Waveguide 

The structure of tapered waveguide is shown in section 3.3, it is simply a rib waveguide 
except it’s width varies firom one end to another end, due to variation in waveguide 
width the propagation constant of rib structure changes while etched section remain 
unchanged. For the calculation of propagation constant, the whole length of waveguide 
is divided into number of thin slices of thickness dy, where y varies from 0 to L. 



Figure 5.17 Electric field (E) distribution in taper waveguide with respect to length (L) 

and width (w). 



Now by the effective index method, effective indices and propagation constants for each 
slice can be calculated, finally it is clear that the effective index and propagation 
constant increases with decreases in waveguide width. The affect of this increase in 
propagation constant and effective index is that the signal weakly guided at smaller 
width side with respect to larger width side, fig.5.17 shows the electric field profile for 
given waveguide structure, at the end A signal is launched to the waveguide when it 
propagates, the value of propagation constant {fi) increases, at the end B propagation 
constant attains maximum value while A attains the minimum value. 



Figure 5.18 Distribution of power (?) in taper waveguide with respect to length (L) 

and width (w). 

The above fig.5.18, shows the power distribution characteristic for tapered waveguide, 
with increase in propagation constant with distance (y) its profile changes at end B 
amplitude of power is very much reduces from the previous value, it means most of the 
power lost out of the rib structure. 



5.6 Couplers 
5.6.1 Parallel Coupler 

In the section 4.1.1 we have discussed the structure of parallel coupler, it is simply two 
parallel waveguides of different width wi and which is separated by a distance ‘r’, 
the geometrical structure is shown in fig.4.2, 



Figure 5.19 Distribution of power (P) in symmetrical parallel optical coupler. 

apply the effective index method to calculate the effective indices of the actual 
waveguide structure, now the real waveguide structure converts into 5 -layers 
waveguide, which is shown in fig.4.3. The fig.5.19 represents the power distribution 
characteristic of parallel coupler, at the end y=0 energy is launched to the waveguide A 
and it will completely coupled to the waveguide B after certain length Lc, this length is 
known as coupling length (Lc). The simulation result shown in fig.5.19, it is clear that 
coupling length Lc for given parallel coupler is 300 pm, which is near to 260 pm 
calculated by Eq. (4.24). 




5.6.2 Tapered Coupler 

5.6.2. 1 Single Tapered Coupler 

Tapered coupler is similar to parallel coupler except its one waveguide is tapered in 
shape, by effective index method the effective indices of actual waveguide can 
calculate. In the earlier section 5.5 we have discussed the method how to calculate the 
effective index of tapered waveguide, similar method we will apply here to calculate the 
effective index of tapered waveguide. 



Figure 5.20 Distribution of power (P) in single tapered optical coupler. 

Now the actual waveguide structure is converted into 5-layer waveguide, shown in fig. 
4.4, which can be easily solved by characteristic matrix method given in section 2.1. 
The simulation result, which is shown in above fig. 5.20 it is clear that coupling length 
Lc for given tapered coupler (0=0.2°) is 260 pm, which is much less then 300 pm for 
parallel coupler, the plot of coupling length (Lc) against tapering angle (0) as shown in 
fig. 5.21, it is clear that effective index of tapered waveguide changes with change in 
length. 





Figure 5.21 Variation of coupling length (Lc) with tapering angle (0) for single tapered 

coupler 

5.6.2,2 Double Tapered Coupler 

In the double tapered coupler, both waveguides are tapered at one end and separated by 
a constant distance ‘ 5 ’, 
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Figure 5.22 Distribution of Power (P) in double tapered coupler 



it any signal is launch in waveguide A then it will coupled in waveguide B after 
coupling length (Lc), tte fig. 5.22 represents the power distribution curve in 
double tapered coupler, the coupling length (Lc) vs tapering angle (0) 
characteristK of double tapered coupler for different separation ( 5 ) of 
waveguides is shown in fig. 5.23. 
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Figure 5.23 Variation of coupling length (Lc) with tapering angle (6) for double tapered 

coupler. 
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Figure 5.24 Schematic representation of multilayer asymmetrical coupler in which, both 

waveguides are tapered. 



The insertion Loss in directional coupler =10 log— dB 

Pq 

Where P;=Input power launched in waveguide A at j=0 
Po=Output power from waveguide B at>=Lc 
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Chapter 6 


Conclusion 

Three- layer semi-infinite dielectric waveguide have been studied earlier, but this 
method is not suitable for N-period stack or multilayer waveguide so special methods 
has been developed for the determination of guided mode and dispersion equation, the 
numerical technique, characteristic matrix method, is introduced here for the 
determination of guided modes in multilayer optical waveguide, in this method 
transverse electric (TE) and transverse magnetic (TM) where deal separately. This 
characteristic matrix method is a generalized method, which is applicable for any kind 
of waveguide structure; by the strict investigation into dispersion relations for 
multilayer systems we can plot the normalization propagation constant {b) against the . 
normalized frequency (V) characteristic, by this characteristic we can calculate the 
effective index and the number of propagated modes in the waveguide structure, once 
we get the effective index of guided structure we can plot the electric and magnetic field 
distribution in the waveguide structure. The application of this method has been carried 
out for the mode field and intensity distribution of planar waveguides. Coupling 
between optical waveguides can be treated as a five-layer medium problem. Such a five- 
layer medium consists of two guiding layered with higher index of refraction, spaced by 
a central layer with lower index of refraction. A MATLAB program has been developed 
to implement the characteristic matrix method. The simulation result represents the 
power distribution characteristic in symmetrical, and asymmetrical (Single and Double 
tapered coupler). Finally it is concluded that coupling length reduces with increase in 
tapering angle (0) for both single and double tapered coupler and for same tapering 
angle (0) the coupling length reduces with decrease in separation (s) of two waveguides. 
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